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Abstract 


In this paper, we study an extension of Riemann—Liouville fractional 
derivative for a class of Riemann integrable functions to Lebesgue measur- 
able and integrable functions. Then we used this extension for the approxi- 
mate solution of a particular fractional partial differential equation (FPDE) 
problems (linear space-time fractional order diffusion problems). To solve 
this problem, we reduce it approximately to a discrete optimization prob- 
lem. Then, by using partition of measurable subsets of the domain of the 
original problem, we obtain some approximating solutions for it which are 
represented with acceptable accuracy. Indeed, by obtaining the suboptimal 
solutions of this optimization problem, we obtain the approximate solutions 
of the original problem. We show the efficiency of our approach by solving 
some numerical examples. 
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1 Introduction 


Fractional calculus is a generalization of classical calculus, which introduces 
derivatives and integrals of fractional order. Major reviews on the concepts 
and history of fractional calculus can be found in the book of [21]. 

Fractional differential equations (FDEs) have been found useful and ap- 
plicable in science and engineering, such as physics, control theory, biology, 
finance, biomechanics, and electrochemical processes (see [2,27] for more 
details). Most of these FDEs do not have analytic solutions; so there are 
some papers dealing with their numerical solutions, for example, predictor- 
corrector method [29], the Adomian decomposition [8], and the variational 
iteration method [19]. Podlubny in [22,23] used the matrix expression to 
unify the formula, where the fractional derivatives are derived from a finite 
difference. Liu also did some interesting works on the numerical approach 
for the fractional differential equations in [2,25] and a method based on 
collocation using spline functions given by [4]. 

On the other side, different models using FPDEs have been proposed 
in [8,10], and there has been significant interest in developing numerical 
schemes for their solution. The great difficulty to obtain the numerical solu- 
tions of such problems is their solutions, which often are represented in terms 
of Mittag—Leffler functions, where these functions are in the form of series 
and their computation is not so easy. One method for solving the FPDE is 
pursued in the recent paper of [LI]. They transformed this partial differen- 
tial equation into a system of ordinary differential equations, which is then 
solved by using backward differentiation formulas. In another very recent pa- 
per, [26] proposed a finite difference method for the fractional diffusion-wave 
equation. [5] considered an implicit numerical scheme for fractional diffu- 
sion equation. The authors of [8] proposed a method for the solution of time 
fractional order partial differential equations by converting it into a nonlin- 
ear programming problem. [16] constructed an efficient spectral method for 
the numerical approximation of the space-time fractional diffusion equations, 
and [5] developed a finite element method for the time fractional Fokker— 
Planck equation. We refer the interested readers to see [2] for more infor- 
mation. 

Based on the above review, we observe that the numerical methods for 
the FPDEs are abundant and when the function is only integrable are very 
limited. As far as we searched the related literature there were no reports on 
similar results as we are investigating in this paper. In this work we focus 
on a novel applicable approach with a better computational cost based on 
optimization problem. We consider a class of fractional convection-diffusion 
equations with measurable (or Lebesgue integrable) functions on the real line 
R and obtain a numerical approximation for the Riemann—Liouville deriva- 
tive of Lebesgue measurable functions. Using the approximated functions 
whenever needed by a joint application of minimization the total error, we 
transform the original FPDE into a discrete optimization problem. By ob- 
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taining the optimal solutions of this problem, we obtain the approximate 
solution of the original problem. The results are more accurate and more 
useful than the ones introduced in [Z| and should be very interesting for 
applications in sciences and engineering. 

The discussion in the rest of paper will be as follows: in the next section, 
we introduce several definitions for different types of fractional derivatives 
of Lebesgue (or Riemann) integrable functions and the notation used in the 
numerical approximation for FPDE. The case with the FPDEs is displayed 
in section 3. Furthermore, in this section, we design an efficient approach 
to approximate the Riemann-—Liouville fractional derivative and use it in our 
numerical approach. Finally, we will give some numerical examples in Section 
4. Conclusions are included in the last section. 


2 Preliminaries 


First of all, let us introduce the properties of measurable functions. 


2.1 Riemann—Liouville fractional extension to Lebesgue 
integral 


Assume that (X, M, j) is a fix measure space [6]. 


Definition 1. If (X,M) is a measurable space and A C X, then the char- 
acteristic function xa of A is defined by: 


ale te A, 
xa(t) = (1) 
0, tZA. 


It is easily checked that 74 is measurable if and only if A € M. A simple 
function on X is a finite linear combination, with complex coefficients, of 
characteristic functions of sets in M. Equivalently, u : X — C is simple 
function if f is measurable and the range f is a finite subset of C. Indeed, 


we have 
n 
U= ) ViX Ai 
i=1 


where A; = u-4(%), 1 <i <n, and range(u)={%1, Y2,---,Yn}- 


Theorem 1. Ifu: X — [0,co) is measurable, there is a monotone increasing 
sequence {un} of simple functions such that 0 < uy < ug < +++ <u, Un > u 
pointwise, and u, > u uniformly on any set on which u is bounded. 
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Proof. See p. 47 of [6]. 


Theorem 2. Let u, : X — [0,00) be an increasing sequence of measurable 
functions on X and for a.ex € X, u(x) = limn+0 Un(x). Then: 


u= lim Un: 
xX noo xX 


Proof. See p. 76 of [24]. 


Theorem 3. If un: X — (0,00) is a finite or infinite sequence of measurable 
functions on X and u(x) = Yo, Un(x) for a.ex € X, then fyu= >>, fy un- 


Proof. See p. 51 of [6]. 


In the special case when the ys is Lebesgue measure on R, the integral, 
which we have developed, is called the Lebesgue integral. At this point it 
is appropriate to study the relation between the Lebesgue integrals and the 
Riemann integrals on R. 


Theorem 4. Let u be a bounded real-valued function on [a, b]. 


(a) If u is Riemann integrable, then u is Lebesgue measurable and fru = 
le d] ue 


(b) wis Riemann integrable if and only if {x € [a,b]: u is discontinuous at x} 
has Lebesgue measure zero. 


Proof. For more details of this theorem, see p. 57 of [6]. 
We shall generally use the notation r u(t)dt for Lebesgue integrals. Now, 
to compute the Lebesgue integral of u, we define the following partition on 


[0, 1]. 


Definition 2. We define the regular partition P,, on [0,1] by: 


Pp = {0 = 60,61,-"° oi), 


where 69 < 01 < +++ < bn. Here, according to the regularization of the 
a 

partitions for each i =0,1,...,n, we have 6; = — and the partition norm is 
n 

defined as: 


\|Prlloo = max {|6; = ial} = : 


In the following, we introduce the definitions of fractional integral and 
derivatives of Lebesgue measurable and integrable functions that can be seen 
as a generalization of the classical derivative. 
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2.2 Fractional calculus 


There are several different ways to define the fractional derivatives, and 
the most commonly used fractional derivatives are the Grunwald—Letnikov 
derivative, the Riemann Liouville derivative, and the Caputo derivative. For 
the definitions of fractional derivatives and some of their applications, see [I4]. 

Let [a,b] € R. We will denote the space of all measurable and Lebesgue 
integrable real functions defined on [a,b] by Ly [a, b]; that is, 


b 
Z,{a,8) = {u: lullasioo =f |u(t)|at < oo}, 


and the space of all measurable and essential bounded real functions defined 
on [a, 6] is denoted by L™{a, }]. 


Definition 3. [£3]. The fractional integral operator of order a > 0 of a 
function u(.) € L1([a, 6], R”), is defined as: 


Amis a / pu (2) 


where T is the classical gamma function and the integral of (2) exists and is 
in the sense of Riemann integral. 


We may claim that the fractional integral (2) is a linear operator on the 
space of Lebesgue measurable functions; that is, if f,g € L1[a,b] and C isa 
constant, then: 

alf (Cf +9) =Calp(f) + al? (9). 


Proposition 1. The fractional integral operator If, 0<a< 1, ts bounded 
or equivalently continuous in Ly[a, |: 


(b= a)° 


alf(u)[lza(a.4) S Pallullzatas), Pa = sa 
| t (u)|lo,5 so] = lull oa 0] T(a+ 1) 


Proof. See its proof in [13]. 


Definition 4. Let u(.) € Ly[a,b].The Riemann-Liouville fractional deriva- 
tive (RLFD) of order a> 0 is defined as: 


Qa d™ m—a _ 1 d™ ‘ m—-a-l 
aDiult) = Fa oll ul) = po ae f I aelar, 8) 


wherem—-1l<a<mandmeN. 


The above fractional derivative is equal to the fractional derivative of 
Riemann-integrable functions just on the compact subsets of R. On the other 
hand, there are some functions that are only measurable and not Riemann 
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integrable but its derivations compute from (B). In this case, consider the 
following function: 


1, t € Qn (0, 1], 


=| (1 
—1, t€ Q°n (0, 1]. 


We know that f(.) is not Riemann integrable function but its Lebesgue in- 
tegrable. Now, we want to compute the fractional derivative of f(.), for 
0<a< 1, as follows: 


i. af a 
Ta) oa (t—7) °f(r)dr 
1 d 


= tarycar— f tor) “dre. 
rd —_ a) = Fowa! ) ae, 


As, m([0, #]) = m([0, 4] M Q] + m([0, t] N Q°) = t, we will have 


oDi f(t) = 


-1 d ‘ 
= ae oT t—1)~dr. 
T(1 _ a) dt eed: ) 


According to Theorem W, the above statement is equivalent to Riemann in- 
tegral and the result of this computation has been shown in Figure [] for 
a = 0.5, 0.8, 0.9, 0.95. 


oDF f(t) = 


GFD 


Figure 1: Fractional derivative of function (@) with different values of a 
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The following theorem help us to apply a fractional integral over a frac- 
tional derivative. 


Theorem 5. Let a> 0. 


(1) aD? aDu(t) = aD?t*u(t), n € N, where N is the set of natural 
numbers. 


(2) The equality oD? alfu(t) = u(t) holds for every u(.) € £1(0,1). 
(3) For u(.) € £1(0,1), n = [a] +1, tf olf °u € AC”1[0, 1], then 


qr-k-1 


n—-1 —a a—k-1 
a ea) =) FSI (SE eli 
k=0 


in which AC|a, b] denotes the space of absolutely continuous real func- 
tions on [a,b], and 


n—1 
AC™|a, b] = {u : [a,b] > R: oo 


AC{a, bl}. 


In particular, if0<a<1 and ,D?u(t) € Lila, bj, then 


(t—a)*-1 


al; aD? u(t) = u(t) on (a) 


( alt *u(a)). 


Proof. See pp. 241-256 of [3]. 


Now, we introduce the definitions of Caputo fractional derivative as fol- 
lows. 


Definition 5. [£3]. Leta > 0, n= [a], andu € AC"[a,b]. The Caputo 
derivative of order a > 0 is defined as 


CDPult) = of Fu) = py f t-re Sar. 


In the space of the functions belonging to AC™[a, b] the following relation 
between Riemann-Liouville and Caputo derivative holds [9]. 


Theorem 6. For u € AC”[a,b], n = [a], the Riemann—Liouville derivative 
of order a of u exists almost everywhere and it can be written as 


n—-1 r—a k-a 
oDeu(t) = Deut) + > pe ula) (5) 
k=0 


In order to simplify our discussion, throughout this paper, we will con- 
sider the Riemann—Liouville definition, since most useful tools have been 
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established by using the Riemann—Liouville definition. It is worthwhile to 
note, by virtue of (4), that for the homogeneous condition considered here 
the Riemann-—Liouville definition coincides with the Caputo version. 


3 Numerical approximation of FPDEs 


Consider the space-time fractional convection-diffusion equation in one di- 
mension as follows: 
O°u(a,t) OFu(a,t) 
Ot OxP 


=f(x,t), 2,t€[0,1],a€(0,1],@e(1,2], (6) 
subject to the following initial and boundary conditions: 


u(z,0) = g(x), x € [0,1], (7) 


u(0,t) = u(1,t) = 0, t € [0,1], (8) 


where g is given a smooth function and f is a measurable function on (0, 1]. 
Here, oe is defined as the RLFD of order a and for all x € [0,1] it is given 
by: 

O° 1 d 


wet) = aay af u(a,s)(t—s)-“ds t€ [0,1], (9) 


and for all t € [0, 1], a is the RLFD of order { that is given by: 

08 t. -a 

— ,t) = ———~ —> t)(x — €)'*d. 0,1). 10 

Fee) = aes | WE N(e- Ode we [0.1 (10) 
Now, we assume that the integrals in equations (Q) and (0) are well-defined. 
Also, suppose that the function u satisfies the Carathéodory conditions with 
respect to L,[0,1], which means, a non-negative measurable function is pos- 
itive on a set of positive measure or, equivalently, the following conditions 
hold: 


(C1) For each x € [0,1], the mapping t > u(a,t) is Riemann integrable. 
(C2) For a.e t € [0,1], the mapping x +> u(a,t) is continuous on (0, 1]. 
(C3) Let u(a, t) : [0, 1] x [0,1] — [0, co) be a bounded function a.e on L;(0, 1]. 


A standard assumption of continuity of the inhomogeneous term can be 
replaced with that of satisfying the Carathéodory conditions with respect to 
L,(0, 1]. 

With specific a and 8 in equation (@), when a = 1 and 6 = 2, we have 
the following classical diffusion equation: 
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Ou(az,t) 0?u(a, t) 
Ot Ox? 


= f(x; t); x,t € [0,1]. (11) 


In fact the time derivative and space derivative of integer order in (I) can 
be obtained by taking the limit a > 1~ in (@) and 6 > 27 in (I). 

In the following, we present a numerical approximation for the Riemann— 
Liouville derivative of Lebesgue measurable functions such that these opera- 
tors simplify to the classical RLFI, and Riemann-—Liouville fractional deriva- 
tives (RLFD) and also the numerical method that gives an approximate so- 
lution to the fractional convection-diffusion equations with measurable (or 
Lebesgue integrable) functions on real line R, studied in this paper. 

In the analysis of the numerical approach that follows, we will assume 
that this space-time fractional convection diffusion equation has a unique 
and sufficiently smooth solution. To establish the numerical approximation, 
and in order to simplify the notation, we consider 


t 
Te" (t) ay u(x, s)(t — 8) “ds, 0<a<l. (12) 
0 
Agta) = f wed@-O Pde, 1<8<2 (18) 
0 
Let t, = nAt, n = 0,1,...,K, where At := fa is the regular partition 
of time interval and x; = jAw, j = 0,1,...,N and Ag := 4 is the regular 


partition of space interval and K and WN are two positive natural numbers. 
To motivate the construction of the scheme, we now define two operators 
that will be focus on them. We use the following approximation at t, for all 
O0<n<k, 


(2 te") (tn) = [TE Gn) — TE" a)], (14) 


and the following approximation at x;, for al0O<j<N, 


(5 AB") (a) = oa [AB s—1) — 2AB(@s) + AR" Cei4s)]- (15) 


For each t,, and x,;, respectively, we need to calculate T27"(t,,) and A aie): 
So, first, we compute (2) by approximating u(a,s) at the fixed constant x 
by simple functions. Therefore, for each s € [0,1], we define U(x,s) as an 
approximation of u(x, s) at the fixed constant x, such that 


U(z, s) _ doula, ti)xas(s); (16) 


in which A; is a sequence of disjoint sets in [0,1] and ¢; = ¢, $= 1, 2, 2g lh. 
Now, we prove the uniformly convergence theorem of our approach. 
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Theorem 7. Let u(x,t) satisfy conditions (C1) — (C3), and let U(ax,t) be 
defined in (1G) for the fixed constant x. If K tends to infinity, then U(a,t) 
tends to u(x,t) uniformly on (0, 1]. 


Proof. The proof is obtained from Theorem [L. 


This theorem means where we choose K as a sufficiently big number, 
closeness U(x, t) to u(x,t), at the fixed constant x, is independent of t € [0, 1]. 
So, by using the approximation (1G) for u(x, s) and suppose that s € [0,t,], 
we will have 


Ta(tn) = i "Yo ule, ti)x4s(6) én — 8) ds, 


that is an approximation to T*™(t,). So, 


Ia(tn) = So ulests) [ta a)*xad(s)de = Youla,t) f (ta—s) "ds. 


i=l 


We denote [0,tn] MA; with A;,,, and note that every bounded Riemann 
integrable function defined on a bounded interval is Lebesgue integrable and 
the two integrals are the same (Theorem M). Then we have 


Ig (t,) = a >» u(x, ti) — (tr - sie = : S- u(a, ti) Wi, 


Ain l-a 


and t; = + for 7 = 1,2,...,n. 

Finally an approximation for (14) is given by 4 [Ia(tn) — Jo(tn—1)]; that 
is, 

(<2) (ty) Pu 1 S- u(a, 4) (Wi = wie) ); 
dt a (1 = a)At a an an-1 

Denote U? as an approximation to the value u(xj;,t,) on grid point (xj, tn). 
Then we take the following approximation for time fractional derivative (Q) 
appeared in problem (@) for the fixed instant «; = jAz € [0, 1], 


ae 4 u(ax;,tagi®, 0<a<l, 
sour = (18) 
u(xj,tn) — u(v;, tr—1) pers 
At : , 


(0) _ wi — wi? ,, and then we have 


an 1, ai,n—1? 


in which q 
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(19) 


in 


ie {" i—1)t-* —2(n—i)-% + (n-i41)-% i<n-l, 


Incidentally, we find that las . ) are bounded for all a € [0,1] and allz > 1, 
as proven in the following Lemma. 


Lemma 1. For all a € (0,1) and alli > 1, it holds Jas? | <c, where c is 
dependent on n and a. 


(x) 


Proof. First, for a = 1 a direct calculation shows that q;,, is a constant, for 


all 1 > 1. Now we prove the lemma for a € (0,1). It can be verified that 
Ain-1 © Ain Cc (0, tr] = m(Ajn-1) < m(Aj,n) < ty. 


In addition, It is well known that if A C B and f is a measurable function 
on A and B, then [, f < J;, f, and by consequence, we have 


O< Wie 1s wi < es 


So the proof is completed. 


Lemma 2. Let U(t) = ~ u(t;)xa,(S), as explained in (IG), be a sum of 
non-negative monotone increasing of simple functions that approximates the 
function u(t) € L,[0,1]. Then 6f(U(t)), which is achieved by (8), tends to 
Dé ,u(t) as K tends to infinity. 


Proof. Using Lemma [], and taking into account the fact that the phrase 
(t —s)~*,0<a< 1, is a continuous and integrable function and y4,(s) isa 
measurable function on [0,¢], and then according to the convergence theorems 
of measurable functions (see Theorem 2 and Theorem of this paper), the 
proof is obvious. 


Similarly, we approximate u(€,t) at the fixed instant t by simple functions 


U(é,t) = DL, ula, t)xa,(8), 


. (20) 
UN, Bi = [0, 1], Bi NB; = 0; i F Js v= WN 


where B,,; = [0,2,;]N Bi, By = =: x], and #7; = Z. Now, we can approxi- 
mate the space fractional derivative (10) for the fixed instant t, = nAt € [0, 1] 
to form 
(Ax)~? 
T(3— 8) 
u(£j-1,tn) — 2u(x;j, tn) + u(tj41, tn) 
Ax? 


Sue, ta, be p>, 


cu = (21) 


B=2, 


12 S. Soradi Zeid, A. V. Kamyad and S. Effati 


in which g{ =W, - 2W) + Wi), and WY) = | - (a; - 5)?-9] 


ij ij-1 


Bij 
and we have 


(j —i + 2)?-8 — 3(7 —¢ +1)? 8 + 3(j — i)?-8 — (f -i-1)2-8, i<j -1, 
gf) = ¢ -3 427-8, t= J, 


1, t= jH+1. 
(22) 
Now, for obtaining the solution of space-time fractional convection-diffusion 


equation (§)-(8), we apply (8) and (21) to (@) and reach the following nu- 
merical method, 


OfUP O20) =F (Gin) (23) 
U9) Ue = 0, fala ale (24) 
So, if u(x,t) is a solution of system (@)-(8), equivalently it is a solution of 


system of equations (23)—(24). Hence, we mention the following main theorem 
of this section. 


Theorem 8. The truncation error for the numerical approximation (B23) - 
(24) is of order O(Ax)?-8 + O(At)!~*. 


Proof. Let u(x,t) be a solution to the fractional partial differential equation 
(@) in which satisfies the conditions (@) and (&). Then we have 


o° 1 d vw 
pa U(e, tn) = Td a) a‘? (tn) 
= rasa ae (Tee = 14) +e1(tn), (25) 


in which €(tp) = O(At). Let us define the error E(t), such that, 
Ta" (tn) — Ta (tn-1) = La(tn) — La (tn-1) + Fr (#). (26) 
So, we have 


Z H(t) tea) + ricw mE + O(AL). 


peat 'n) = Tay Re 


Now, by using Lemma [I] we have E(t) = O(At)?~®. So it follows that: 


o* 1 1 


pe tltste) = Fay (Ja(tn) = Ta(tn-1)) a. O(At)!-* + O(At). 


T(1—a) 
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Therefore, 


OX 


0 a n —a 
(Sat) itn) = PUP + O(A8)-*. (28) 


Similarly, for the fractional derivative of order 8, we have 


la - = 
(S25) (a;, tn) = 62UP + O(Aa)?-?. (29) 
Finally, by using the numerical method (23)—(24), we will have 


UF — uF loo = O(At)“* + O(Ax)?-*, (30) 


and the proof is completed. 


In the remainder of this section, it will be shown that equations (23)—(24) 
can be converted to an optimization problem. For this suggested approach 
we need the following theorem. 


Proposition 2. If h : [a,b] > [0,00) is a measurable function, the necessary 
and sufficient condition for ft h=0 is that hh =0 a.e., on [a,b]. 


Proof. See p. 51 of [6]. 


So, according to Proposition 2, necessary and sufficient condition for 
u(x,t) to be a solution of system (23))—(24) is that the optimal solution of the 
following problem is zero (see Theorem 1 of [I] for convergence): 


1 pl 
min ‘| | 
o Jo 


By trapezoidal method and using the ending point in any subinterval for 
approximating integrals and using approximations (TX) and (21) for (BI), we 
obtain the following discretized problem: 


6°U (a, t) — 68U (x, t) — f(a, t)|dtde. (31) 


‘ihe n jt 
min 7 S- S- Jao ule; tal _ by u(ai,tn)al? — cf (xj, tn)|, 
j=ln=1 i=l i=1 


(32) 


where a = (Ax)T(3 — 6), b = (At)°T(2 — a), c = ab, and a ae are as 
before. Now, we convert the problem (B62) to a linear programming problem 
with the following change of variables (see Theorem 3 and Lemma 2 of [I] 


and [28] for more details): 
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K 
> (7 +67) (33) 


ln=1 


AtAaz 


Cc 


N 
min 
Uu A 


J 


n j+1 
s.t. ay) ui - by ung? =6f @jjtn) 17 — OF 
i=1 i=l 


Finally, by obtaining the solution of problem (83), we recognize the value of 
unknown variables. 


4 Numerical examples 


In this section, we give some numerical examples and apply the presented 
approximation for solving them. These test problems demonstrate the valid- 
ity and efficiency of this techniques. 


Example 1. We compute 9D@x(t), with a = 0.5, for x(t) = t* by ap- 
proximation (8). 
The exact formulas of the derivatives are derived from 
D9 (t8) = Ti(s+1) ps 0.5 
: T'(s+1-—0.5) ; 
Figure J shows the results. Since the exact solution for this problem is known, 
we compute the approximation error by using the maximum norm for each 


Kk. Assume that Z(t,), n = 1,...,K, are the approximated values on the 
discrete time horizon t,,...,t«%. Then the error is given by 
E = max (|2(tn) — 2(tn)|). (34) 
n 


In the case of a = 0.5, with various choices of kK, the maximum absolute 
errors are computed by equation (84) and shown in Table IL. 


Example 2. Consider the following fractional differential problem: 
oD? a(t) = g(x(t), t), (35) 


where 0 < a < 1, 2(0) = 0 and g(z,t) = oe -e 


a(t) +t? -t. 
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Analytic 


(42 n 4 4 | 
0.45 0.455 0.46 0.465 0.47 0.476 0.48 


Figure 2: Analytic solution and numerical approximation ((Q), with various choices of 
K, for Example 1 


Table 1: Maximum absolute error for Example 1 


Kk Eappro 

5 2.65757 x 10-8 
10 6.47115 x 10-4 
20 4.09595 x 1074 
40 9.77932 x 10-5 
50 2.27352 x 1075 
100 5.04375 x 107 


The exact solution of this equation is x(t) = t? — t. We can rewrite the 
equation (BH) as follows: 


oD? x(t) — g(x(t),t) = 0. (36) 
Also, we have from equation (80) that 

oD x(t) — g(x(t), t)| = 0. (37) 
By using again Proposition , we conclude that if x(t) is a solution of equation 


(B12) with the initial condition «(0) = 0, equivalently it is a solution of the 
following optimization problem: 


ath i | oDea(t) — 9(e(t),t) at (38) 
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Then by discretization the integral as before an using equation (19) for ap- 
proximate of 9Dfa(t), we simplify obtained problem (8) as follows: 


N K k 
; At) 0 
min i= 5 2i(tp) ae (39) 
“j=l k=1 p=1 
242-% 
— (ate(S# p20 a) (wi(te) — + te))], 


where t, and re are the same as before. Now, for solving the minimum 


problem (BY), we can change it to a linear programming problem like (83). In 
Figures Bl and @] we compare the exact solution with numerical approximation 
(BY) for two different values of kK, N = 100, and a = 0.5. 


a1 T T T T T T T T T 
* approximation 
exact 


Figure 3: Analytic solution and numerical approximation (GY) for Example 2 for K = 40 


Table 2 shows the exact solution and the approximate solution for equa- 
tion (B54) by using problem (89) for K = N = 100 and a = 0.5,0.99. The 
results compare well with those obtained in [20]. From numerical results 
we can indicate that the solution of FDE approaches to the solution of in- 
teger order differential equation, whenever a approaches to its integer value. 


Example 3. We consider the space-time fractional diffusion equation (@) 
with the following initial and boundary conditions: 


u(x, 0) = 0, 0<a<1, 


u(0,t)=u(1,t)=0, te [0,1], 


where 
t?-sin(2rxz) + 4n*t? sin(2r2). 
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0.1 T T T T T T T T T 
* approximation 
exact 


Figure 4: Analytic solution and numerical approximation (9) for Example 2 for K = 100 


Table 2: Numerical values for Example 2 when a = 0.5, 0.99 


t Lappro(a _ 0.5) Lexact Lappro(Q = 0.99) Lexact 
0.0 0.000000 0.000000 0.000000 0.000000 
0.1 -0.089978 -0.090000 -0.089586 -0.090000 
0.2 -0.159889 -0.160000 -0.159688 -0.160000 
0.3 -0.209891 -0.210000 -0.209707 -0.210000 
0.4 -0.239974 -0.240000 -0.239787 -0.240000 
0.5 -0.249896 -0.250000 -0.249738 -0.250000 
0.6 -0.239998 -0.240000 -0.239795 -0.240000 
0.7 -0.199879 -0.210000 -0.209830 -0.210000 
0.8 -0.160109 -0.160000 -0.159897 -0.160000 
0.9 -0.096390 -0.090000 -0.100098 -0.090000 


The exact solution of this problem when 8 = 2 is given by u(az,t) = 
t?sin(2rx). The spatial and temporal meshes are taken uniform; that is, 
At = %, Av = =. From the Table B and Figures 6, 6, and MZ, we get that the 
approximate analytical solution is consistent with the exact solution when 
a—>1, 6 > 2 and different values for K and N. Therefore, the validity of 
our numerical methods is confirmed. 


Now, consider the vector U(Ax) = (Up,U1,--: ,Un), where U; is the 
approximate solution, for 2; = jAx, 7 =0,1,...,N at a certain time t, and 
u(Ax) = (u(xo,t),...,u(an,t)), where u is the exact solution. The absolute 
error is defined as follows: 


Ag) — (Ae )llea = mane [elegy t) UG, (40) 
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Table 3: Maximum absolute error (@M) for the numerical scheme of Example 3 at t = 1 
for different values of a and 8 = 2 with various choices of At and Ax 


At Ax Eappro 
1/5 1/5 0.00847 
a=05 1/10 1/10 0.00431 


1/15 1/15 5.32982 x 10-4 
1/20 1/20 7.05965 x 10-5 
1/5 $1/5 0.00480 
a=09 1/10 1/10 1.64153 x 10-4 
1/15 1/15 7.19164 x 1075 
1/20 1/20 7.39623 x 10~® 
Ws - ifs 0.00536 
a=0.95 1/10 1/10 1.00857 x 10-4 
1/15 1/15 3.38148 x 1075 
1/20 1/20 2.41435 x 10-6 
is 1/5 0.00568 
a=0.98 1/10 1/10 2.31643 x 10-4 
1/15 1/15 4.62727 « 10-* 
1/20 1/20 5.70822 x 10-6 
1/5 1/5 0.00579 
a=0.99 1/10 1/10 1.14488 x 10-4 
1/15 1/15 9.19805 x 1075 
1/20 1/20 2.82570 x 10-6 


Figure 5: The evolution of u(a,t) for anomalous diffusion coefficients a = 0.5, 8 = 2 for 
Example 3 
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Exact Solution 


Figure 6: The evolution of u(z,t) for anomalous diffusion coefficients a = 0.9, 6 = 2 for 
Example 3 


Exact Solution 


Figure 7: The evolution of u(a,t) for anomalous diffusion coefficients a = 0.98, B = 2 
for Example 3 


Table @ presents a comparison between the absolute error of our method and 
the method presented by Neville et al. (2011) in [@] with different space and 
time steps. It should be noted that our numerical results, quickly converging 
to the exact solution with a lower division, are more accurate than with the 
results obtained in [6]. 
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Table 4: Comparision of absolute error for the numerical scheme of Example 3 at a = 0.5 
and § = 2 with various choices of kK, N 


kK ON our method Kk WN © estimated in [7 
5 5 0.000360 4 100 0.000926 
10 10 0.000133 8 100 0.000337 
15 15 5.12841x10-® 16 100 0.000120 
20 20 8.33727x 10-6 32 100 4.184925 x 107° 


Example 4. we consider the following space-time fractional diffusion prob- 
lem: 


d°u(a,t) (2.8) u(x,t) _ 0.8 Ou 
ate 2 axb = q(z,t) —2 Ox’ (x,t) € (0, 1] x (0, 1] 
(41) 
with the initial condition: 
u(z,0) = #7(1— 2), x € [0,1], (42) 
and the boundary condition: 
u(0,t) = u(1,t) = 0, t>0, (43) 
where 
2n?(1 — x)? 


The exact solution of to this FPDE for a = 0.8 and 6 = 1.5 is given 
by u(z,t) = 27(1 — x)(1+ #7). Table B shows the absolute errors at time 
t = 1, between the analytical solution and the numerical solution obtained 
by applying the method discussed in this paper and the method presented 
in [BOj. It should be noted that our results are very closely identical with 
other results. 


Table 5: Comparision of absolute error for the numerical scheme of Example 4 for a = 0.8 
and 8 = 1.5 at t = 1 with various choices of K and N 


Kk  N- our method — estimated in [BO] 
10 10 1.23 1073 1.26 x 1073 
20 20 5.86 x 10-4 6.74 x 1074 
40 40 2.87x 1074 3.48 x 1074 
80 80 83x 10-5 8.6 x 1078 
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5 Conclusion 


FDEs and FPDEs have caught considerable attentions due to their various 
applications. Since there is no systematic method to derive the exact solu- 
tion of these equations, in this paper we partly solved a class of FPDE and 
obtained some approximating solution for it by provide an efficient approach 
based on measurable functions. For this purpose, at first, we present an 
approximation for the fractional derivatives of measurable functions. Then, 
by using this approximation and application of minimization the total error, 
we transform the original FPDE into a discrete optimization problem such 
that the optimal solutions of this problem is the approximate solution of the 
original problem. The suggested method represents a unifying approach for 
the solution of partial differential equations of fractional order. 

From the numerical examples and the results that are compared with the 
exact solutions and with the other methods, it is shown that, as the number of 
discrete points, K and N, in proposed approach of this paper was increased, 
the solutions converged to the exact solutions. In addition, as the value of 
qa approaches one and ( approaches two, the numerical solutions for both 
the FDE and FPDE approach the analytical solutions for a = 1 and 6 = 2. 
Since the proposed approximation of this paper is based on the minimization 
of total error, it is clear from the results that there is no difference between 
exact and approximate solution in point to point case. 
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